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We calculate I C (T) and R„ for both coherent and incoherent electron tunneling across a c-axis 
break junction between two v — s, d x 2_ y 2-wa,ve layered superconducting half spaces, each with c- 
axis bandwidth 2J. Coherent quasiparticle tunneling only occurs for voltages V < 2J/e, leading 
to difficulties in measuring R n for underdoped samples. The coherent part of I c {0) is independent 
of A 1/ (0) for J/A„(0) << 1, and can be large. Our results are discussed with regard to recent 
experiments. 
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It is presently possible to prepare high quality o 
axis Josephson junctions of Bi 2 Sr2CaCu208+<5 (BSCCO). 

Early thin film atomic layer-by-layer molecular 
beam epitaxy (ALL-MBE) preparations of trilayer junc- 
tions of BSCCO separated by a thin layer of Dy-doped 
BSCCO were found to have values of the product I c R n 
of the critical current times the normal resistance at tem- 
perature T = consistently about 0.5 mV, with I c and 
R n varying greatly. E| Recently, they prepared a sin- 
gle Josephson junction within a single unit cell of un- 
derdoped BSCCO, and reported I c R n « 5 - 10 mV. @] 
Also, a blunt point contact tip pressed onto a BSCCO 
surface often resulted in an apparent c-axis break junc- 
tion. Overdoped junctions typically had I c R n ~ 2.4 
mV, well below the Ambegaokar-Baratoff (AB) result. 
However, two underdoped break junction samples had 
I c R n (0) w 15 — 25 mV, apparently in excess of the ex- 
pected BCS-derived AB result. Furthermore, exceed- 
ingly clean c-axis break junctions were prepared by cleav- 
age and subsequent refusion of BSCCO, with or without 
a twist about the c-axis. || Although I c (T)R n was not 
presented in these works, we anticipate that such mea- 
surements will be available shortly. 

Here we consider tunneling across an untwisted c-axis 
break (or single intrinsic Josephson) junction which is 
much less conductive than the bulk, intrinsic junctions 
between neighboring Cu02 layer pairs. We find that for 
underdoped samples, standard measurements of R n at 
voltages V > 2A(0)/e, where A(T) is the superconduct- 
ing order parameter (OP) amplitude, can be unreliable, 
since they do not measure the coherent processes, which 
can dominate at V = 0. Hence, the apparent large val- 
ues of IcRn reported for underdoped samples could be 
questionable, flp) 

We assume a c-axis break junction between two un- 
twisted half-spaces of cross-sectional area A, each con- 
sisting of A >> 1 identical clean superconducting lay- 
ers separated a distance s apart. We label the up- 
per (it) and lower (£) half-spaces by fj, = u, £, and in- 



dex the layers in each half-space by j = 1, . . . , N, with 
j = 1 being the layer in each half-space adjacent to 
the break junction. We allow ^ MiJj(T (k) [^ J CT (k)] to 
annihilate [create] a quasiparticle with spin a = ±1 
and two-dimensional (2D) wavevector k on the j th layer 
within the /i th layered half-space. Within each layer in 
the /x th half-space, the quasiparticles propagate with en- 
ergy dispersions ^(k) = e^oC 4 ) — Ep relative to the 
Fermi energy Ep [for free particles, e^o(k) = k 2 /(2m M )], 
and interact with intralayer BCS-like pairing interac- 
tion A M (k, k') = X^voPviMfvifa'), where v = s,d, 
^s(</>k) = 1 and (pd{4>k) = v / 2cos(20k) are the eigen- 
functions for s and d x 2_ y 2-wsLve intralayer pairing, re- 
spectively. We only consider here the purely s-wave and 
rf-wave cases X ns0 ^ X n > d0 = and X nd0 ^ X^ s0 = 0. Be- 
tween neighboring layers in the fi th half-space, the quasi- 
particles tunnel with matrix element J^/2. The c-axis 
resistivity p c (T) above the transition temperature T c sug- 
gests the limits J^/T c » 1 and J n /T c « 1 apply to 
overdoped and underdoped materials, respectively. S 

In addition, we take the single particle tunneling 
Hamiltonian Ht across the break junction to be 



Ht = E 7k,k^Ii,,(k)V^i, CT (k') + H.c. 



(1) 



k,k' 



which transfers a quasiparticle from the j = 1 layer in 
the I half-space to the j = 1 layer in the u half-space, 
and vice versa; 7k,k' = T£, f We- set h = c = Ub = 1. 

We assume both coherent and incoherent break junc- 
tion tunneling. The spatially constant coherent tunnel- 
ing preserves the intralayer wavevectors, k = k'. How- 
ever, pure spatially random incoherent tunneling assumes 
k and k' are independent of each other, || which allows 
no d-wave incoherent Josephson tunneling. Hence, to al- 
low for a finite (albeit extremely small) amount of d-wave 
incoherent Josephson tunneling, we assume to second or- 
der inT k , k ,, §§ 



7k,k'?k',k' 



ASkj 



|T | 2 A5 k , k , +/ i nc(k-k / ) , (2) 



1 



where 



/i„ c (k - k') 



1 



2nN 2D (0) 



1 2cos[2(0 k - 



(3) 



1/rxd « l/r±„ and N 2D (0) = [N iDu {0)N 2 Dt(P)] l/2 is 
the geometric mean 2D density of states; for free parti- 
cles, N 2 d(0) = w/(2tt), where m = (m^mf) 1 / 2 . In Eq. 
(g), < • ■ • > denotes a 2D spatial average. 

For intralayer pairing in a half-space with one conduct- 
ing layer per unit cell, the regular and anomalous temper- 
ature Green's functions G M (k, u>) and F^j (k, u>) , where 
u> represents the Matsubara frequencies, for propagation 
within layer j in the fi half-space, explicitly depend 
upon j. [p|JlO[| However, the OP A p is independent of 
j. §] Nevertheless, A M (c/> k ) = A jUI/ (T)yi I/ (0 k ) implicitly 
depends upon v = s,d. jll| 

For purely s-wave incoherent tunneling between 3D su- 
perconductors, I c R n is independent of the properties of 
the junction jj]. However, in our model these quanti- 
ties must be evaluated separately. R n is found from the 
quasiparticle current I qp to leading order in , |fl2f 



4e 
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k,k' 



k,k' 



de[/(e»)-/( Q )]x 



*[G u ,i(k ) -ie u )]3[G M (k , ) -ie<)], 



(4) 



^° h (^) = W° h (^) = 64G/(3^) 7 6(1 - 7 2 )Q( 7 ), (7) 



where C = 2e\T \ 2 N 2D (Q), Q( 7 ) = (l + | 7 |)[(l+ 7 2 )i;(fc)- 
2\ 7 \K(k)}, k = (1-| 7 |)/(1 + | 7 |), 7 = eV/(2J), 
K{z), E(z) are standard complete elliptic integrals, and 
0(z) is the Heaviside step function. In the bulk state 
approximation, Q( 7 ) is replaced by (3/8)-ftT(fc)/(l + | 7 |), 
which leads to a spurious, non-Ohmic InV* dependence 
of l/iC h as7->0. 

Although Igp h (V) is Ohmic for V — > 0, it is non-Ohmic 
for finite V, and vanishes for \eV\ > 2J. Thus, for 
\eV\ > 2J, the only quasiparticle tunneling across the 
break junction process allowed is incoherent. This result 
arises from the geometry: the half-spaces are assumed 
to be layered, each with c-axis bandwidth 2J. For 3D 
systems, there is no such limitation upon Iqp h - 

Thus, the quasiparticle current consists of two parts, 



Igp(V) = I™ + I c q ; h (V) = V/R n {V). In the inset 



of Fig. 1, we have plotted R /R n (V), where Rq = 
J/[87r 2 e 2 7V 2 D(0)|7^| 2 ], as a function of eV/2J. Note that 
one requires the break junction conductance to be much 
less that the conductance across neighboring layers in 
each half space. This implies that both 1/t± s and \%\ 2 / J 
are small with respect to J 2 t\\, where 1/t|| << J is the 
small intralayer scattering rate. However, this does 
not restrict the relative magnitudes of Ig° h and P n nc 



qp 



where f{x) is the Fermi function, e u = e, ei = e + eV", and 
the G Mj i(k, — iEp) are obtained from G AIj i(k, oj) by the 
analytic continuations to — > — ie^. Since the tunneling 
takes place between the j = 1 in the two half-spaces, the 
only relevant wavevectors are 2D. Hence, we set — > 

AN 2D (0) d^o Jo * #k/(27r). '] 

We consider separately the coherent and incoher- 
ent processes, and separately the G Ali i(k, to) as eval- 
uated exactly for the layered half-spaces, and as ap- 
proximated using the bulk layered states. In the bulk- 
space treatment, we assume G Mi i(k, w) « G Al ,b(k, lo) 
= Jo ^/[iuj-^o - J^cosz], or G ft6 (k,u/) = l/R p (iu>), 
where R^{z) = [{z — jUio) 2 — J 2 ] 1 ^ 2 depends upon k 
only through £ M o(k). pf However, when one properly 
takes account of the surface at the weak break junc- 
tion, H G^i = ^ sin 2 z/[iu - ^ - J^cosz], or 
G M ,i(k, to) = E^iJ), where 



E ll (z) = 2/[z-^ + R tl {z)]. 



(5) 



Using either the bulk or half-space states and the identity 
S-oo - /( e + eV )} = eV > the incoherent quasipar- 

ticle c-axis break junction tunneling current is Ohmic, 



rinc 

qp 



V/RZ = 2e 2 VN 2D (0)/r ±s 



(6) 



For the coherent quasiparticle c-axis break junction 
tunneling current, we only consider tunneling between 
identical materials with Ji = J u = J, etc., and obtain 
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FIG. 1. / c c ^(J,0)//^ h (0,0) and P c %{ J, 0)// c c °^( J, 0) are 
plotted versus log 10 [J/A„(0)] for v — s,d. Inset: Plot of 
Ro/Rn(V), where R = J/[87r 2 e 2 7V 2 o(0)|To| 2 ], and its bulk 
space approximation versus eV/2J. The curves are shifted by 
J/[47r 2 |T | 2 T Xs ] values of 0.05, 0.1, 0.15, and 0.2 for clarity. 

In the superconducting state, I c for Josephson tunnel- 
ing across the break junction between arbitrary layered 
half spaces is given to lowest order in 7k t k' by j?]] 

j c( T ) = ^rE r k ,k^M(k,^;>>), (8) 



Lj,k,k' 



where quite generally = -A M (0 k )9r M /.D M , = 
[exp (ik+s) - exp (iks)] /(iJ^) and = [|A M | 2 + 

uj 2 ] 1 / 2 . j|] The quantities exp(ik±s) are obtained from 



2 



the equation J M cos(fc±s) = — £^o(k) ± iD^, M which 
leads to T M = S M (iD M ), where is given by Eq. (|). 

In the bulk state approximation for F^ i, is replaced 
by R~ 1 (iD f i). Note that as J M — > 0, the half-space and 
bulk expressions both reduce to the familiar 2D form. 

The incoherent part of the break junction I c between 
arbitrary layered superconductors with v = s, d is 



2eN 2D {0)Tir 



En 



2ji 



2irD„ 



(9) 



Note that A /X (<^k) implicitly depends upon v. Equation 
(Q) is obtained using either the bulk or half-space states. 
For s-wave pairing, I™(T)R™ C equals the AB result. |] 
For A„ = A e = A, j£„ c (0) = TreiV^^A^O)^/^, 
where A s — 1, and = 0.153. For this case, 
4 n s c (T)/4 n s c (0) and 4 n d c (T)/4 n d c (0) are plotted in Figs. 
2 and 3, respectively. ' We note that I™l{T)/P™(T) cx 

TLs/Tld << 1. 

For coherent c-axis break junction Josephson tunneling 
between identical materials, we drop the subscripts /i, 
noting that A(</> k ) = A„(T)yv((/> k ), etc. Writing T in F x 
in integral form as above Eq. (||), and performing two 
integrals analytically, we have 



ACT 



2tt 



* dzX(z) 

o 



(10) 



| l/ 2 y/ \ -4 

A(z) = sin z 



where W(z) = [l + J 2 sin 2 z/D 2 ] 
sin 2 z/[l + W(z% and C = 2e|T | 2 ^ 2D (0). Using the 
bulk states, I%£ b (J,T) is obtained from Eq. @ by re- 
placing X(z) by sin 2 z/8. In the limit T — > 0, we set 
«(0 k ) = \J/A„(0)MM\, /?(0k) = 1(1 + [1 + « 2 ] 1/2 ), 
and Icf(J,0) = CY„[J/ A v (0)] reduces to 



/™V,o) = 



4C 



1 + 8/3 
24/3 2 



ln/3 

™2 



(11) 



For the bulk states, /£„ b (J, 0) is obtained from 
i™( J, Q) by replacing the quantity in large brack- 
ets by (8a) _1 sinh -1 a. In Fig. 1, we have plot- 
ted 7 c c ° h (J,0)// c c ° h (0,0) and P c %( J, 0)//™ h b (0, 0) versus 
log 10 [,// Aj,(0)] for s- or d-wave OPs. The most surpris- 
ing point is that for small J/A y (0), 1^(0) = 1, so that 
I™v(J,ty -> C, independent of A„(0). 

From Fig. 1, 7^ h (J, 0) is slightly more sensitive to J 
than is I^° h (J, 0). In addition, for intermediate values of 
J/A„(0), the respective s- and d-wave bulk curves closely 
approximate the half-space curves obtained by reducing 
J/A v {0) by the constant factor l/s/2. For J/A„(0) >> 
1, however, there are distinct differences between the 
bulk and half-space curves. Whereas the correct — 
CY u [J/A v (0)] -> 16CA V (0)B V /(3J), where B s = 1, and 
Bd = 2\/2/vr, the approximate 7c° h h (</, 0) spuriously re- 
duces to CB v [A v {Q)/J\\j^J/D v K v {Q)}, where D s = 1, 
D d = y/2/e. 



It is interesting to compare the coherent and inco- 
herent results for identical half-spaces. At T = 0, the 
V = T q ; h (0)/I™ cx \T \ 2 t ±s /J, whereas for T > T c , 
J/T c is the relevant quantity for the half spaces. j(| 
Since A„(0) « T c , for J/T c « 1, / c c ° h ( J, 0)/7^ c (0 ) cx 
\%\ 2 t±„/T c . For J/T c » 1, r c f(J,0) and 7^(0) both 
cx A„(0), but 7 c c ° h (J,0)/4 n , c (0) cx \%\ 2 t^/J. These re- 
sults lead to the curious conclusions that for J/T c << 1, 
the underdoped normal state tunneling is incoherent, the 
T = quasiparticle break junction tunneling could be ei- 
ther coherent or incoherent, but d-wave break junction 
pair tunneling would be mainly coherent. On the other 
hand, for J/T c >> 1, the overdoped normal state half- 
space tunneling would be coherent, but the T — quasi- 
particle break junction tunneling and the s-wave break 
junction pair tunneling could be incoherent! 

In the limits J M — > 0, one can evaluate the J = 0, 
T = limit of the coherent part of I C {T) from Eq. 
(§) as a function of r = A u /A e , obtaining I c coh (0,0) = 
2CVln(r)/(r 2 — 1), where C is given following Eq. (0). 
In the limit r — ► 1, 7™ h (0, 0) — > C. For either two s-wave 
or two <i-wave superconductors, r is independent of </> k . 

In Fig. 2, we have plotted I c co s h ( J, T)// c co s h ( J, 0), as a 
function of T/T c , for tunneling between two layered s- 
wave superconductors. Typical curves with J/A s (0) = 
0,2,10 are shown, along with the AB curve, Eq. (J|). 
For J/A s (0) = 100, the curve is almost identical to the 
AB curve. Using the bulk states doesn't change these 
curves very much, except for large J/A s (0). Clearly, 
I C , S {J,T) 1 7 CiS (J,0) is rather indistinguishable from that 
of AB, independent of the microscopic details. 




FIG. 2. Plots of J c c ° h (J,T)/7 c c ° h (J,0) for J/A»(0) = 0,2, 10 
and of 4" C (T)/4" S C (0) (AB) versus T/T c , for tunneling be- 
tween identical s-wave half-space superconductors. 

In Fig. 3, we have plotted 7™ h ( J, T)//™ h ( J, 0), as 
a function of T/T c . Typical curves with J/A(0) = 
0,0.5,2,10 are shown along with the ci-wave analog of 
AB, 4 n d c (T)/^ n d c (0). Note that the magnitude of I££(T) 
is very small, due to the factor of 1/rxd- Unlike the s- 
wave curves in Fig. 2, the I c °}{ J, T)//™ h ( J, 0) curves 



3 



with small values of J/A^O) are distinctly linear at low 
T, and are thus distinguishable from the AB curve, in 
Fig. 2. However, for J/A d (0) » 1, J c c ° h ( J, T)/I c c f( J, 0) 
and I™d(T)/I™j(0) are are nearly indistinguishable from 
the AB curve. 
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FIG. 3. Plots of I c °i{J,T)/r°^{J,0) for J/A d (0) = 0,0.5, 
2, 5, 100 and of 7^(T)/7^(0) versus T/T c , for tunneling be- 
tween identical d-wave half-space superconductors. 

Summarizing our results, we have for v = s,d, 

where Z v = A v t± s /t± l , and ry = 32|7o| 2 tj_ s /(37t J). Since 
I c is measured at V = 0, one requires R n (0). For over- 
doped samples, J/T c >> 1, and I CAJ (J,0)R n (0) reduces 
to ttA^O)^ + r]B„]/ {2e[l + ??]}. For both v = s,d, 
this is proportional to A„(0), and nearly independent of 
the break junction properties for s-wave superconductors. 
For c?-wave superconductors, one requires a substantial 
/"J 1 in order to obtain a non-negligible I c R n - However, 
for underdoped samples, J/T c << I, and the situation is 
far more complicated. First of all, one cannot measure 
R n (0) in the usual way, since the coherent contribution, 
which can be large at V = 0, vanishes for eV/2A„(0) > 1. 
Moreover, / C J/ (J, 0) is dominated by coherent tunneling 
and independent of A„(0) for |7o| 2 >> irA v A„(0)/2T_ Ll/ , 
which is especially likely for c?-wave superconductors. 

In conclusion, we found that for c-axis break junc- 
tion tunneling between two layered superconductors, a 
crossover from incoherent quasiparticle to coherent pair 
tunneling can occur. This greatly complicates the de- 
termination of R n (0) , the coherent part of which cannot 
be seen from measurements with eV/2A„(0) > 1, un- 
less J >> A v (0), which corresponds to overdoped sam- 
ples. For underdoped samples, I c R n values tend to be 
overestimated. The approximate bulk electronic states 
lead to correct incoherent, but incorrect coherent, tun- 
neling results. Incoherent <i-wave pair tunneling leads 
only to very small I c R n values. For coherent pair tun- 
neling, both s-wave and d-wave pair tunneling are large 



in magnitude for small J/A ly (0), and cross over to the 
AB form for large J/A„(0). The T-dependence of coher- 
ent d-wave tunneling is distinctly different from that for 
AB for small J/A„(0). Thus, accurate measurements of 
the T-dependence and magnitude of I c R n in such break 
junctions could give important information regarding the 
questions of the order parameter symmetry and of the 
coherence of the pair tunneling. 
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